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Abstract. 

Magnetic field relaxation is determined by both the field's geometry and its topology. For 
relaxation processes, however, it turns out that its topology is a much more stringent constraint. 
As quantifier for the topology we use magnetic helicity and test whether it is a stronger condition 
than the linking of field lines. Further, we search for evidence of other topological invariants, 
which give rise to further restrictions in the field's relaxation. We find that magnetic helicity is 
the sole determinant in most cases. Nevertheless, we see evidence for restrictions not captured 
through magnetic helicity. 
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1. Introduction 

Geometry and topology of magnetic field lines fundamentally affect their dynamics 
(IWoltierlll958l: I Arnold Il974l : iRuzmaikin et a/.lll994l : iTavlorl 11974 : iDel Sordo e7aZll2010l : 
Yeates et aZ.ll2Q10l . 1201 lh 7 For instance, strongly tied field lines give rise to strong current 
sheets which then facilitate magnetic reconnection under which field lines brake and 
connect in a different way. Reconnection for its part, can give rise to ejections of plasma, 
which is of particular interest in the case of our Sun. 

While the field's geometry has often been appreciated, its topology has received less 
attention. Loosely speaking, topology determines the field's linkage, while geometry its 
configuration in space. Any two field configurations which are topologically different 
cannot be transformed one into the other without breaking field lines, i.e. reconnection. 



2. Magnetic helicity 

Magnetic helicity density is the scalar product of the magnetic vector potential A and 
the magnetic field B, i.e. 

h = A-B. (2.1) 
Its integral over a closed or periodic system, the total magnetic helicity, 



H 



I 



A-B dV, 



(2.2) 



is a conserved quantity in ideal MHD and in the limit of vanishing magnetic resistivity 
( WoltierlfT958h . 

Topologically speakin g, H is a quant ifier for the mutual linkage of magnetic flux tubes 
and their internal twist (Moffatt Il969h - Twisted and linked fields are severely restricted 
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Figure 1. Iso surfaces of the magnetic energy for the initial magnetic field configurations. 
Arrows denote the direction of the field. The left configuration is non- helical, while the center 
one is helical. The right configuration was used as control run and is not helical. 



in their dynamics, in particular their relaxation. Arnold! ( 1974 ) first quantified this re- 
striction in the realizability condition 



E(k) =2\H(k)\/k, 



(2.3) 



with the spectral magnetic energy E(k), the spectral magnetic helicity H(k), and the 
wave number k. It gives a lower bound for the magnetic energy in the presence of magnetic 
helicity. In its picturesque interpretation as linking of flux tubes it becomes clear why 
magnetic helicity imposes restrictions on the magnetic field decay given by equation 
(j2.3p . During relaxation, mutually linked field lines cannot freely evolve without magnetic 
reconnection. As long as reconnection is not aided by strong inflows of magnetic fields 
into the reconnection zone, it will not occur fast enough for any appreciable field change 
or energy loss. There exist, however, field topologies of linked magnetic field lines which 
are not helical and for which equation (j2.3p has no effect. 

In a first work we investigate whether the field's topology, as it is given by the linking 
and twisting of field lines, is the dete rmining factor in rela xation, or whether the magnetic 
helicity content is the key quantity ( Del Sordo et a/.ll2010l ). From the plethora of possible 
magnetic field configurations one of the simplest examples is chosen, which is a triple 
ring configuration of interlinked flux tubes (Fig. [I]). 

We solve the resistive MHD equations for a viscous, compressible and isothermal gas 



d 

— A = U x B-rjvoJ, 
at 

D 9 1 

—U = -CgVlnp+ - J xB + F v 

S-hip= - V-E7, 
Dt 



(2.4) 
(2.5) 
(2.6) 



with the velocity Z7, the molecular magnetic resistivity 77, the susceptibility in vacuum 
/io, the electric current density J = V x B/po, the isothermal speed of sound c s , the 
density p and the advective time derivative D/Dt = d/dt + U ■ V. Viscous effects are 
caught in F v [ sc = p~ x V • 2^pS, where v is the kinematic viscosity, and S is the traceless 
rate of strain tensor with components S^- = \{u%^ + u j,i) ~ ' U . Commas denote 
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Figure 2. Magnetic field lines for the two interlocked triple ring configurations after 4 Alfvenic 
times of resistive decay. The non-helical initial configuration (left panel) loses its shape quicker 
than the helical configuration (right panel). 
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Figure 3. Normalized magnetic energy evolution for the helical linked rings (solid line), 
non- helical linked rings (dashed line) and unlinked rings (dotted line). 



partial derivatives. Initial magnetic fields represent either of the three configurations in 
Fig. [TJ while the initial velocity vanishes in the whole domain and the initial density is 
unity. Boundary conditions are chosen as periodic in order to conserve magnetic helicity. 

From the time evolution of the magnetic field lines it becomes clear that linking alone 
cannot hinder the fast decay of the magnetic energy (Fig. [2j left panel). In the presence 
of magnetic helicity, however, the decay is slowed down considerably (Fig. [2j right panel). 
Our control setup with non-interlinked flux tubes shows an energy decay characteristics 
which is very close to the interlinked non- helical field (Fig. [3]). The helical configuration, 
on the other hand, shows a much slower decay rate. From this we conclude that magnetic 
helicity, rather than linking of flux tubes, determines the field's dynamics. 

In the same fashion we test the importance of the field's linkage and knottedness for 
non-helical configurations with other highly non-trivial topologies. Initial fields are the 
IUCAA knot (Fig. [4j left panel), which, in the Alexander-Briggs notation, is the 8_18 
knot, and the Borromean rings (Fig.HJ right panel), after the emblem of the north Italian 
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Figure 4. Iso surfaces of the initial magnetic energy for the IUCAA knot (left panel) and the 

Borromean rings (right panel). 




Figure 5. Comparison of the magnetic energy evolution for different initial field 

configurations. 



aristocratic house of Borromeo. We compare the magnetic energy evolution with the triple 
ring configurations and find that both, the IUCAA knot and the Borromean rings, show 
an intermittent power law in the energy decay (Fig.EJ). This allows for speculations about 
higher order topological invariants, which are non-zero for tho se field configurat i ons an d 
impose additional restrictions on the field's dynamics; see e.g. iRuzmaikin et~al\ ( 1994 ). 



3. Beyond magnetic helicity 

There exists an infinite number of topological invariants for three-dimensional vector 
fields. Applications on MHD have been, nevertheless, very limited and successful in only 
a few attempts. Two invariants of third and forth order in B are finite for the Borromean 
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rings ( Ruzmaikin et a/1 ll994). Their definition is, nevertheless, limited to distinct flux 
tubes which do not overlap. In resistive MHD magnetic field diffuses and an initially 
confined field will occupy the whole space. 

One way around this hitch is by using the fixed point index (see e.g. iFrankell ( 2004f ): 
Yeates et all (|2010h ). which in turn, is only applicable to fields with a preferential direc- 
tion, like in toroidal fields, or a field with a positive z-component. In the latter case one 
can define a mapping between the bottom and top boundaries by tracing the field lines, 
resulting in the field line mapping, 



-> K 2 , 
(x,y) F(x,y), 



(3.1) 
(3.2) 



with the initial point (x,y) at z = 0. Note that F(x,y) is bijective. 

Fixed points are those points for which F(x,y) = (x,y). There can be infinitely or 
finitely many, with at least one fixed point. Considering the fixed point's neighborhood 
in the xy plane we can determine its sign. Depending if F x (x, y) > x and F y (x, y) > y, a 
color is assigned; for F x (x, y) < x and F y (x, y) > y a different color is assigned, likewise 
for the other cases. The sequence of these colors around the fixed point determines its 
sign U, where U £ { — 1,1} . Summing over all fixed points yields the fixed point index 



(|Brownlll97il iFrankefeOOl 



(3.3) 

which is a conserved quantity in ideal MHD (|Brownlll97lh . 

Simulations using the fixed point index as constraining quantity were performed by 
Yeates et all ( 2010[|2Q11 ). where they observed a constraint relaxation of magneti c fields. 
Their equilibrium state turned out to be of higher energy than that proposed by iTavlorl 
(|l974h . 



4. Conclusions 

From resistive MHD simulations of relaxing interlinked magnetic fields it becomes ap- 
parent that magnetic helicity, rather than actual linkage, determines the field's relaxation 
properties. The decaying IUCAA knot and Borromean rings show some intermittent de- 
cline speed for the magnetic energy. This suggests that there might be higher order 
topological invariants, which impose restrictions on the field's dynamics. An example of 
such invariants is the fixed point index, which is conserved in ideal MHD and is shown 
to impose further restrictions on relaxation. 
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